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ABSTRACT 

V 393 Scorpii is a bright Galactic Double Periodic Variable showing a long photometric cycle 
of « 253 days. We present new VIJK photometric time series for V 393 Scorpii along with 
the analysis of AS AS V-band photometry. We disentangled all light curves into the orbital and 
long cycle components. The ASAS V^-band orbital light curve was modeled with two stellar 
components plus a circumprimary optically thick disc assuming a semidetached configuration. 
We present the results of this calculation, giving physical parameters for the stars and the 
disc, along with general system dimensions. Our results are in close agreement with those 
previously found by Mennickent et al. (2010) from IR spectroscopy and the modeling of the 
spectral energy distribution. The stability of the orbital light curve suggests that the stellar 
+ disc configuration remains stable during the long cycle. Therefore, the long cycle should 
be produced by an additional variable and not-eclipsed emitting structure. We discuss the 
evolutionary stage of the system finding the best match with one of the evolutionary models 
of van Rensbergen et al. (2008). According to these models, the system is found to be after an 
episode of fast mass exchange that transferred 4 from the donor to the gainer in a period 
of 400.000 years. We argue that a significant fraction of this mass has not been accreted by 
the gainer but remains in an optically thick massive (~ 2 M©) disc-like surrounding pseudo- 
photosphere whose luminosity is not driven by viscosity but probably by reprocessed stellar 
radiation. Finally, we provide the result of our search for Galactic Double Periodic Variables 
and briefly discuss the outliers j3 Lyr and RX Cas. 

Key words: stars: early-type, stars: evolution, stars: mass-loss, stars: emission-line, stars: 
variables-others 



1 INTRODUCTION 

V393 Scorpii is one of the Galactic Double Periodic Variables 
(DPVs), a group of interacting binaries showing a long photomet- 
ric cycle lasting roughly 33 times the orbital period (Mennickent 
et al. 2003, Mennickent & Kolaczkowski 2009, Michalska et al. 
2009, Poleski et al. 2010). DPVs have been interpreted as semi- 
detached interacting binaries with ongoing cyclic episodes of mass 
loss into the interstellar medium (Mennickent et al. 2008, Mennick- 
ent & Kolaczkowski 2010). The 253-d long photometric cycle of 
V 393 Scorpii was discovered by Pilecki & Szczygiel (2007) after 
inspection of the ASAS catalogue for eclipsing binaries with ad- 
ditional variability. The lUE-UV properties of V 393 Scorpii were 
studied by Peters (2001) who found evidence for a hot temperature 
region produced by the tangential impact of the gas stream into the 
gainer photosphere. This region should be the origin of the superi- 
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onized lines observed in the UV, like N V, CIV and Si IV, that are 
likely produced by resonance scattering in a plasma of temperature 
T ~ 10^ K and electron density ~ 10^ cm"^ (Peters & Polidan 
1984). The star was also studied by means of multi-epoch high- 
resolution IR spectroscopy by Mennickent et al. (2010, hereafter 
MIO), who also studied broad-band photometry and lUE ultravi- 
olet spectra. After summarizing the available literature of this ob- 
ject, these authors argued for a semidetached B3 + FO binary with 
masses 8 M0 and 2 M© for the gainer and donor (hereafter also 
called primary and secondary, respectively) and orbital separation 
of 35 -R0. Most remarkably, MIO found evidence for large mass 
loss through the Lagrangian L3 point during epochs of long cycle 
minimum and claim that their observations suggest that the mass 
loss producing the long cycle is probably concentrated in equato- 
rial regions. 

In this paper we refine stellar and system parameters of V 393 
Scorpii by fitting the light curve with a multicomponent model in- 
cluding a stationary circumprimary disc. A detailed report of the 
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Table 1. Summaiy of photometric data analyzed in this paper. ASAS data 
can be found at http://www.astrouw.edu.pl/asas/ 



Observatory 


Telescope 


Filters 


N 


HJD-range (-2450000) 


La Silla 


REM 60 cm 


V 


359 


4498.8391-5372.6688 


La Silla 


REM 60 cm 


I 


340 


4498.8402-5391.7978 


La Silla 


REM 60 cm 


J 


210 


4498.8393-5362.8939 


La Silla 


REM 60 cm 


K 


169 


4502.8541-5369.8908 


Las Campanas 


ASAS-array 


V 


903 


1950.8740-4973.8121 



observations used in this paper is given in Section 2, our results 
are presented in Section 3, a detailed discussion of these results is 
given in Section 4 and our conclusions are presented in Section 5. 




orbital phase 



2 OBSERVATIONS 

We obtained V,I,J,K band images with the robotic 60-cm REM 
telescope operated at La Silla, Chile, during three seasons in 2008- 
2010. Bias and flat calibration images were regularly obtained, 
and de-biasing and flat-fielding for all science images were done 
with standard photometric IRAlQ tasks. Instrumental magnitudes 
were calculated using aperture photometry and also by applying the 
differential image analysis technique (DIA, e.g. AUard & Lupton 
1998). Differential magnitudes were calculated between the target 
and nearby stars used as comparison stars. The DIA technique pro- 
vided the less noisy light curves and are shown in this paper. We 
also analyzed V-band magnitudes available in the public ASAS 
archiv^ Typical one-sigma error for ASAS and REM photometry 
is 0.03 and 0.05 mag, respectively. This is the kind of statistical 
error provided in the rest of the paper for some modeled or fitting 
quantities. A summary of all photometric observations analyzed in 
this paper is given in Table 1 . 
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3 RESULTS 

3.1 Light curve disentangling 

In this paper we use the ephemeris for the main orbital pho- 
tometric minimum provided by Kreiner (2004), namely To — 
2 452 507.7800 + 7.7125772_E. We separated the light curve of 
V393 Sco into a long and a short period component. For that we 
used an algorithm especially designed to disentangle multiperiodic 
light curves through the analysis of their Fourier component ampli- 
tudes. The method works as follows. A main frequency /i (usually 
the orbital one) is found with a period searching algorithm applied 
to the light curve. A least square fit is applied to the light curve con- 
sidering a fitting function consisting of a sum of sine functions of 
variable amplitude and phases representing the main frequency and 
their additional significant harmonics. After this first fit, the resid- 
uals are inspected for a new periodicity /2 (the long periodicity in 
DPVs). This new periodicity and their harmonics are included in 
the new fitting procedure. Finally, we obtain the best representation 
for the light curve based on a sum of Fourier components of fre- 
quency /i, /2 and their significant harmonics. Data residuals with 

^ IRAF is distributed by the National Optical Astronomy Observatories, 
which are operated by the Association of Universities for Research in 
Astronomy, Inc., under cooperative agreement with the National Science 
Foundation. 
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Figure 1. Orbital light curve (up) and long-cycle light curve (down) of V393 
Scorpii at three bandpasses. Ephemeiis are given in the text. In the lower 
light curve circles represent REM-/f and pluses ASAS-V data. 

respect to the second and first theoretical light curves form the or- 
bital and long cycle photometric series, respectively. 

This deconvolution method was applied to the light curves of 
V 393 Sco at the observed bandpasses yielding a unique solution 
consisting of the orbital light curve plus an additional smooth 
variability with period 253 days. A phase dispersion minimization 
analysis (Stellingwerf 1978) of combined ASAS-V and REM-V 
data yielded for the long-cycle the following ephemeris for the 
light curve maximum: 

HJDmaa:.long = 2452522(±1) + 253(±4) * E (1) 

Inspection of the orbital and long light curves at different 
bandpasses reveals that (Fig. 1): (i) the rounded shape of the or- 
bital light curves at quadratures are consistent with a close binary 
showing proximity effects, (ii) the system shows smooth long-cycle 
variability with larger amplitude in redder bandpasses. The ampli- 
tudes of this long variability in the V, I, J and K bands, derived 
from sinus fits to the REM data are 0.177 ± 0.003, 0.240 ± 0.004, 
0.231 ± 0.003 and 0.263 ± 0.003 mag., respectively, and (iii) V 
and K light curves look different, the maximum of the latter seems 
to be delayed by A^i ~ 0.05 with respect to the former. 

The orbital light curve was examined with the program Pe- 
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riodO'^J We calculated the error in the main frequency of the fit to 
the orbital light curve. The error consistently given by Monte Carlo 
simulations and the method of least squares is 4 X 10"^ Hz. This 
means that the period could drift at most by 4.76 x 10~^ days in 
3023 days (the ASAS-V^ dataset time baseline). This implies that 
a constant orbital period change if present, should be less than 0.5 
seconds per year. 



3.2 Model for a circumprimary optically thick disc 

Here we give a brief description of the disc model that we apply to 
V393SCO. 

The basic elements of the binary system model with a plane- 
parallel disc and the corresponding light-curve synthesis procedure 
are described in detail by Djurasevic (1992, 1996). 

We assume that the disc in V 393 Sco is optically and geomet- 
rically thick. The disc edge is approximated by a cylindrical sur- 
face. In the current version of the model (Djurasevic et al. 2010), 
the thickness of the disc can change linearly with radial distance, 
allowing the disc to take a conical shape (convex, concave or plane- 
parallel). The geometrical properties of the disc are determined by 
its radius {Rd), its thickness at the edge (de) and the thickness at 
the center (dc). 

The cylindrical edge of the disc is characterized by its 
temperature, Td, and the conical surface of the disc by a radial 
temperature profile obtained by modifying the temperature distri- 
bution proposed by Zola (1991): 

T{r)^Td + {Tn-Td)[l-{i^^)T^ (2) 

We assume that the disc is in physical and thermal contact with 
the gainer, so the inner radius and temperature of the disc are equal 
to the temperature and radius of the star (Rh, Th)- The temperature 
of the disc at the edge {Td) and the temperature exponent (ay ), 
as well as the radii of the star {Rh) and of the disc {Rd) are free 
parameters, determined by solving the inverse problem. 

The model of the system is refined by introducing active re- 
gions on the edge of the disc. The active regions have higher local 
temperatures so their inclusion results in a non-uniform distribu- 
tion of radiation. The model includes two such active regions: a hot 
spot (hs) and a bright spot (bs). These regions are characterized by 
their temperatures Tha,bs angular dimensions 6ha,bs and longitudes 
Xhs,b3 ■ These parameters are also determined by solving the inverse 
problem. 

Due to the infall of an intensive gas-stream, the disc surface 
in the region of the hot spot becomes deformed as the material ac- 
cumulates at the point of impact, producing a local deviation of ra- 
diation from the uniform azimuthal distribution. In the model, this 
deviation is described by the angle 6rad between the line perpen- 
dicular to the local disc edge surface and the direction of the hot 
spot maximum radiation in the orbital plane. 

The second spot in the model, i.e. the bright spot, approxi- 
mates the spiral structure of an accretion disc, predicted by hydro- 
dynamical calculations (Heemskerk 1994). The tidal force exerted 
by the donor star causes a spiral shock, producing one or two ex- 
tended spiral arms in the outer part of the disc. The bright spot can 
also be interpreted as a region where the disc significantly deviates 
from the circular shape. 
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3.3 Results of the V-band light curve fitting 

The light-curve fitting was performed using the inverse-problem 
solving method (Djurasevic 1992) based on the simplex algorithm, 
and the model of a binary system with a disc described in the previ- 
ous section. To obtain reliable estimates of the system parameters, 
a good practice is to restrict the number of free parameters by fixing 
some of them to values obtained from independent sources. Thus 
we fixed the spectroscopic mass ratio to g = 0.25 and the donor 
temperature to T2 = 7900 K, based on our spectroscopic study 
(Mennickent et al. in preparation). In addition, we set the gravity 
darkening coefficient and the albedo of the gainer and the donor to 
Ph^c = 0.25 and Ah,c = 1.0 in accordance with von Zeipel's law 
for radiative shells and complete re-radiation (Von Zeipel 1924). 
The limb-darkening for the components was calculated in the way 
described by Djurasevic et al. (2010). 

We treated the rotation of the donor as synchronous (/c = 1.0) 
since it is assumed that it has filled its Roche lobe (i.e. the filling 
factor of the donor was set to Fc = 1.0). In the case of the gainer, 
however, the accreted material from the disc is expected to transfer 
enough angular momentum to increase the rate of the gainer up to 
the critical velocity as soon as even a small fraction of the mass has 
been transferred (Packet 1981, de Mink, Pols & Glebbeek 2007). 
This means that the gainer fills its corresponding non-synchronous 
Roche lobe for the star rotating in the critical regime, and its dimen- 
sions and the amount of rotational distortion are uniquely deter- 
mined by the factor of non-synchronous rotation, which is the ratio 
between the actual and the Keplerian angular velocity. For V 393 
Scorpii we assumed critical rotation for the gainer, i.e. fh = 14.6 
(Model A) but also calculated a model using synchronous rotation 
(Model B) to estimate the effect of the gainer velocity in the system 
parameters. 

The results of the light-curve analysis based on the described 
model of V 393 Sco given in Table 2 basically confirm our solution 
for the stellar parameters given in MIO, based on IR spectroscopy 
and the modeling of the spectral energy distribution. Furthermore, 
now we obtain a more realistic Rh value. Our results obtained con- 
sidering synchronous rotation for the gainer do not differ much 
from the critical rotational case. The fit and stellar and disc di- 
mensions are illustrated in Fig. 2. We note that residuals show no 
dependence on orbital or long-cycle phases. 

We find that the best fit model of V 393 Sco contains an opti- 
cally and geometrically thick disc around the hotter, more massive 
gainer star. With a radius of Rd ~ 9.7Rq, the disc is more than 
twice as large as the central star {Rh ~ 4.4i?Q). The disc has a 
moderately convex shape, with central thickness of dc ~ 2.1_R0 
and the thickness at the edge of de « I.SRq. The temperature of 
the disc increases from Td = 8600 K at its edge, to Th = 16600 
K at the inner radius, where it is in thermal and physical contact 
with the gainer. The temperature profile exponent ar in Eq. 1 is 4.5, 
this means that the effective temperature of the disc is significantly 
higher than the temperature at its edge. 

We were able to model the asymmetry of the light curve very 
precisely by incorporating two regions of enhanced radiation on 
the disc: the hot spot (hs), and the bright spot (bs). The hot spot 
is situated at longitude \hs ~ 320 degree, roughly between the 
components of the system, at the place where the gas stream falls 
onto the disc. The longitude A is measured clockwise (as viewed 
from the direction of the -l-Z-axis, which is orthogonal to the or- 
bital plane) with respect to the line connecting the star centers (-I-X- 
axis), in the range 0-360 degrees. The temperature of the hot spot 
is approximately 20% higher than the disc edge temperature, i.e. 
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~ 10300 K. The hot spot can be interpreted as a rough ap- 
proximation of the hot line which forms at the edge of the gas 
stream between the components (Bisikalo et al. 2003). Although 
including the hot spot region into the model significantly improves 
the fit, it cannot explain the light- curve asymmetry completely. By 
introducing one additional bright spot, larger than the hot spot and 
located on the disc edge at Xbs w 160 degree, the fit becomes much 
better. 



4 DISCUSSION 
4.1 Disc formation 

In a semidetached system with the donor transferring matter onto 
the gainer through an accretion stream, the stream has a distance 
of closest approach rmin from the center of the gainer, given 
approximately by: 



: 0.0488g" 



(3) 



(Lubow & Shu 1975). Using q = 0.25 we obtain rmin = 0.093a, 
comparable to the gainer radius (0.125 ± 0.046 o). The fact that our 
light curve model shows a relatively large disc around the gainer in- 
dicates that: (i) the stream in principle hits the gainer that rapidly 
reaches critical rotation (Packet 1981), it cannot accrete additional 
mass and a disc is built around the star or (ii) the gainer radius is 
smaller than rmin and a disc is formed naturally. At present, uncer- 
tainties in q and Rh do not allow to select between these mutually 
exclusive alternatives. 



4.2 On tlie cause for the long-cycle 

It is notable that the model containing the circumprimary disc fits 
very well the orbital light curve through the whole long cycle. This 
fact suggests that the optically thick disc does not participate in 
the long cycle, at least that part of the disc responsible for the 
orbital photometric variabiUty. We conclude that the long cycle 
is produced by an additional variable and non-eclipsed emitting 
structure, a fourth light in the system. In addition, the fact that 
this fourth light is redder at long maximum and bluer at minimum 
places strong constrains on the possible cause for the long cycle 
variability. 

In principle, it is possible that the long minimum is due to 
obscuration of the system by ejected material through the equa- 
torial plane, as suggested by MlO's observations. This possibil- 
ity is not supported by our observations, since it should produce 
a white (opaque material) or red (gas producing reddening in the 
line of sight and free-free emission outside the line of sight) min- 
imum. On the contrary, our observations indicate that the system 
when brighter is redder, a fact consistent with free-free emission in 
non-obscured ejected circumbinary material. If this material is not 
ejected in the orbital plane at long maximum, then it is ejected at 
higher latitudes. In this sense the inferred variable emitting struc- 
ture is reminiscent of the jets found in /3 Lyr (Harmanec et al. 1996, 
Harmanec 2002.) 

This interpretation conflicts with the MIO suggestion that 
mass loss occurs in the equatorial plane. That position was based 
on the non- variability of UV spectral features through the long cy- 
cle and the evidence for mass flows through the L3 point. However, 
only few UV spectra were available for analysis, and it is possible 
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Figure 2. Observed (LCO) and synthetic (LCC) light-curves of V393 Sco 
obtained by analyzing photometric observations; final O-C residuals be- 
tween the observed and optimum synthetic light curves; fluxes of donor, 
gainer and of the disc, normalized to the donor flux at phase 0.25; the views 
of the optimal model at orbital phases 0.25, 0.50 and 0.75, obtained with 
parameters estimated by the light curve analysis. 
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Table 2. Results of the analysis of V393 Sco V-fllter light-curve obtained by solving the inverse problem for the Roche model with a disc around the more- 
massive (hotter) component in critical rotation regime (Model A) and synchronous rotational regime (Model B). 



Quantity 


Model A 


Quantity 


Model A 


Quantity 


Model B 


Quantity 


Model B 


n 


2268 


MJMq] 


7.8 ±0.2 


n 


2268 


M^[Me] 


7.8 ± 0.5 


E(0 - C)2 


0.5638 


M,[Me] 


2.0 ±0.2 


E(0 - C)2 


0.5621 


MJMq] 


2.0 ± 0.2 


^rms 


0.0157 




4.4 ±0.2 


""rms 


0.0157 




4.1 ± 0.2 




80.0 ±0.2 




9.4 ±0.3 


in 


79.9 ± 0.2 




9.4 ± 0.3 


Fd 


0.55 ± 0.04 


log 9h 


4.0 ±0.1 


Fd 


0.51 ± 0.04 


logflh 


4.1 ±0.1 


Td[K] 


8600 ± 600 


log 5c 


2.8 ±0.1 


Td[K] 


8300 ± 600 


logfle 


2.8 ±0.1 


dcllorb] 


0.04 ± 0.01 


, 

bol 


-3.0 ±0.2 


do [a-oib] 


0.04 ± 0.01 




-2.9 ±0.3 


In, 1 1 


0.06 ± 0.01 


, 

bol 

aorb[R-0] 


-1.4 ± 0.1 


dr. fo, V,l 


0.08 ± 0.01 


, 

bol 
«orb [R0] 


—1.4 ± 0.1 


ax 


4.5 ±0.5 


35.1 ± 0.3 


ax 


4.4 ±0.5 


35.1 ± 0.5 


fh 


14.6 ±0.8 


7^d[R©] 


9.7 ± 0.3 


fh 


1.00 


7^d[R0] 


9.1 ± 0.5 


Fh 


1.00 


de[R0] 


1.3 ±0.3 


Fh 


0.25 ± 0.01 


do[R0] 


1.4 ± 0.3 


Th[K] 


16600 ± 500 


dc[R0l 


2.1 ±0.4 


Th[K] 


16800 ± 500 


dc[R0] 


2.8 ±0.5 


Ahs = Ths/Td 


1.2 ±0.1 






Aha = Ths/Td 


1.3 ±0.1 








19.0 ±2.0 








18.5 ± 2.0 






Ah.n 


324.0 ± 6.0 






Ahs[°] 


325.0 ± 6.0 






''rad[°] 


-5.0 ±5.0 






Srad[°] 


-5.0 ±5.0 






Abs = Tbs/Td 


1.2 ±0.1 






Abs = Tbs/Td 


1.3 ±0.1 






db.l°] 


33.0 ± 6.0 








29.0 ± 6.0 






Ab.n 


162.0 ±9.0 






Abs[°] 


154.0 ±9.0 








9.90 ±0.03 








8.83 ± 0.03 






fie 


2.35 ±0.02 






Oo 


2.35 ± 0.02 







FIXED PARAMETERS: q = Mc/Mh = 0.25 - mass ratio of the components, Tc = 7900K - temperature of the less-massive (cooler) donor, Fc = 1.0 - 
filUng factor for the critical Roche lobe of the donor, fc = 1 .00 - non-synchronous rotation coefficients of the donor, /9h,c = 0.25 - gravity-darkening 

coefficients of the components, Ah.c = 1.0 - albedo coefficients of the components. 

Note: n - number of observations, S(0 — C)'^ - final sum of squares of residuals between observed (LCO) and synthetic (LCC) light -curves, cr mis - 
root-mean-square of the residuals, i - orbit inclination (in arc degrees), = R^/Ryc - disc dimension factor (the ratio of the disc radius to the critical 
Roche lobe radius along y-axis), Td - disc-edge temperature, de, dc, - disc thicknesses (at the edge and at the center of the disc, respectively) in the units of 
the distance between the components, ax - disc temperature distribution coefficient, /h - non-synchronous rotation coefficient of the more massive gainer 
(ratio between the rotational angular velocity and the synchronous rotational angular velocity), Fh = Rh/Rzc - filling factor for the critical Roche lobe of 
the hotter, more-massive gainer (ratio of the stellar polar radius to the critical non-synchronous Roche lobe radius along z-axis for a star in critical rotation 
regime, equal to 1 for critical rotation regime), Th - temperature of the gainer, Ahs.bs = Ths,bs/Td - hot and bright spots' temperature coefficients, Shs,bB 

and Aijs,bs - spots' angular dimensions and longitudes (in arc degrees), 6i.ad - angle between the line perpendicular to the local disc edge surface and the 
direction of the hot-spot maximum radiation, Hb^c - dimensionless surface potentials of the hotter gainer and cooler donor, Ai^^ ^ [Mq], T^j, ^ [R0] - stellar 
masses and mean radii of stars in solar units, log g^^ ^ - logarithm (base 10) of the system components effective gravity, M^^j - absolute stellar bolometric 
magnitudes, Oorb [Rq]. '^d[R©]. depe], dc[R©] - orbital semi-major axis, disc radius and disc thicknesses at its edge and center, respectively, given in 

solar units. 



that the putative variability was hidden by additional orbital vari- 
ability, or alternatively, the long cycle variability disappears at UV 
wavelengths. Much more well-sampled UV spectra are needed to 
illuminate this point. In addition, the inferred mass loss through L3 
around secondary eclipse is still compatible with vertical outflows. 
Furthermore, UV lines do show evidence for such high latitude out- 
flows (MIO). 



constant (e.g. Warner 1995). From the above expressions and using 
the parameters given in Table 2 and T = 8600 K (the temperature 
in the disc outer edge) we obtain M = 3.74 x 10~® Mq yr~^. 

Let's use the formula given by Smak (1989) for the geomet- 
rical thickness of an a-disc at the density level of 10~^° g cm~^: 

^^Q.07[M^sil-i^)°r (6) 



4.3 Mass transfer and mass loss 

In this section we calculate the mass transfer rate (M) for Model 
A using several approaches. The radial temperature structure of a 
steady state accretion disc is given by: 

r(r)=To[4]-^/^[l-[^]V2]i/^ (4) 

where 

rp _ r 3GMhM n/4 ^r-, 

-^0- I 8^0- J{3 J P) 

fl 

where G is the gravitational constant and cr the Stephan-Boltzman 



where Mis is the mass accretion rate in units of 10^® g s ^ . Using 
the average between de and dc as the value for H, we foimd M = 
7.97 X 10"^ Mq yr-i. 

Now we assume that at $0 = 0.25 the system magnitude 
is V= 7.609 (Table 4). From the light curve model the disc flux 
contribution at this phase is 10%. Considering the reddening 
E{B -V)= 0.13 and a distance of 523 pc (MIO) we obtain V^*'" 
= 9.68 and My"" = 1.09. Using the bolometric correction BC = 
-1.50 for T = 17000 K, log g = 3.0, and Z = Zq (Lanz & Hubeny 
2007) we obtain M^^i" = -0.41 corresponding to 115 Lq. If this 
luminosity is powered by accretion, then: 
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Table 3. The parameters of the van Rensbergen et al. (2008) model that best 
fits the V 393 Sco data. The hydrogen and hehum core mass fractions (Xc 
and Yc) are given for the cool and hot star 



quantity 


value 


quantity 


value 


age 


7.00E7 yr 


period 


7.713 d 


Mc 


2.11 Mq 


Mh 


7.49 M0 


Mc 


-9.47E-9 Mg/yr 


Mh 


9.47E-9 Mg/yr 


logTc 


3.92 K 


logTh 


4.32 ft: 


log Lc 


2.61 Lq 


logijl 


3.39 Lq 


Rc 


9.55 R0 


Rh 


3.79 Rq 




0.05 


^ch 


0.63 


Ycc 


0.93 




0.35 



2Bh 



(V) 



(e.g. Warner 1995). Using the above equation we derive M = 4.15 
X 10"^ Mq yr-^ 

If the DPV long cycles are due to recurrent episodes of sys- 
temic mass loss as suggested by Mennickent et al. (2008), and as- 
suming no mass accumulation around the gainer in the long term, 
then AAd = M Piong are ejected from the system every long cycle. 
For the M values given above this means that AAf between 5.4 
X 10~^ and 2.9 x lO"*" Mq are ejected in every cycle. We find 
that in an accretion powered system, the amount of mass lost in ev- 
ery long cycle is not minor and should affect the system evolution 
considerably. 

If the system loses mass due to a spherically symmetric wind 
that does not interact with the companion, then the orbital period 
changes by (e.g. Hilditch 2001): 



Pg 



Mh+Mg 



(8) 



Using our stellar masses and estimates for the mass loss rate 
given above we find changes in the orbital period of 1 to 10 seconds 
per year, easily detectable with the current astronomical instrumen- 
tation. The fact that this variability is not observed argues against 
an accretion disc whose luminosity is driven by viscosity but in 
favor of an extended photosphere radiating by reprocessed stellar 
radiation. This view is also supported by other evidence presented 
later in the next section. 



4.4 Evolutionary stage of V 393 Scorpii 

The comparison of the stellar parameters with predictions of evo- 
lutionary models for single stars with solar metallicity indicates: 
(i) the donor is significantly inflated appearing as an over-luminous 
object for its mass and (ii) the gainer is under-luminous and under- 
heated for the given mass (Fig. 3). The inflated cooler star is ex- 
pected in a semidetached interacting binary with an evolved donor 
whereas the under-luminous hotter star can reveal the presence of 
a massive circumprimary disc. In fact, the loci occupied by the 7.8 
Mq gainer in the HR diagram corresponds to a slightly evolved 
6 Mq star. The figure Rh = 4.4 Rq is also compatible with a 6 
Mq star. This could be possible if a 2 Mq optically-thick disc sur- 



rounds the gainer. In other words, in this view, our dynamical mass 
is overestimated by the presence of a massive discQ 

The average disc electron number density can be estimated 
from: 



(9) 



Using the values given in Table 2 and Mmsc = 2 Mq we derive We 
= 3.2 X 10^^ cm~"^. This density is much larger than that found in 
a normal B-type photosphere and comparable with densities found 
in the stellar interior. 

An alternative explanation involving a less massive disc is that 
the gainer light is in some way obscured by the presence of the 
high latitude wind detected in the UV producing a smaller effective 
temperature and luminosity. Whereas the hypothesis of a massive 
disc explains the small Rh , the obscuration effect does not. 

The hypothesis of a self-gravitating massive disc (0.5 Mq) for 
the DPV-related system P Lyr was considered by Wilson & Terrell 
(1992) but critiqued by Hubeny et al. (1994) who pointed out that 
a massive-disc model requires unrealistically low viscosity, i.e. a 
large Reynolds number, and is likely to be dynamically unstable. 
All these authors assumed accretion-powered discs in their calcu- 
lations, leaving the door open for a structure with a different mech- 
anism for energy generation. Further theoretical work is necessary 
to explore alternative physical scenarios for massive discs around 
normal stars more deeply. 

Now we consider predictions of binary evolution models 
including epochs of systemic mass loss. We inspected the 561 
conservative and non-conservative evolutionary tracks by van 
Rensbergen et al. (2008) available at the Center de Donnees 
Stellaires (CDS) looking for the best match for the system pa- 
rameters found for V 393 Sco. Models with strong and weak tidal 
interaction were studied, although only the latter ones should allow 
critical rotation of the gainer. A multi-parametric fit was made 
with the synthetic (Si.j.k) and observed (Ok) stellar parameters 
mass, temperature, luminosity and radii, and the orbital period, 
where i (from 1 to 561) indicates the synthetic model, j the time tj 
and k (from 1 to 9) the stellar or orbital parameter. Non-adjusted 
parameters were mass loss rate, Roche lobe radii, chemical 
composition, fraction of accreted mass lost by the system and age. 
For every synthetic model i we calculated the quantity xl,j at 
every tj defined by: 



(l/iV)EfeWfc[(S,,,,fc-Ofc)M] 



(10) 



where A'^ is the normalization factor and Wk the statistical weight 
of the parameter Ok, calculated as: 



Wk = v/Ofc/e(Ofc) 



(11) 



where e{Ok) is the error associated to the observable Ok- The 
model with the minimum corresponds to the model with the best 
evolutionary history of V 393 Sco. The absolute minimum Xmin 
identifies the age of the system along with the theoretical stellar and 
orbital parameters. The high accuracy of the orbital period domi- 
nates the search for the best solution in a single evolutionary track. 



* The DPV AUMon suffers of the same problem. Considering Mi = 7 
Mq , log Th = 4.20 and log = 3. 17 found by Djurasevic et al. (2010) the 



gainer turns to be under-luminous. 
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but the others parameters play a role when comparing tracks corre- 
sponding to different initial stellar masses. 

We find the absolute minimum in the weak interaction 
model with initial masses of 6 and 3.6 M© and initial orbital pe- 
riod of 3 days. This model starts Roche Lobe overflow at about the 
same time that core hydrogen burning ends, so it is a borderline 
Case A/B (thanks to Nicki Mennekens for this insight). The aver- 
ages of their parameters are shown in Table 3. 

The corresponding evolutionary tracks for the primary and 
secondary stars are shown in Fig. 4, along with the position for the 
best model for V 393 Sco. We observe a relatively good match for 
the donor star parameters but a mismatch with the gainer tempera- 
ture and luminosity. The higher stellar temperature and luminosity 
indicated by the model could indicate that in practice the gainer has 
not accreted all the transferred mass but part of the mass has been 
accumulated in a massive optically thick surrounding disc. Conse- 
quently, lower temperatures and luminosities are observed for the 
gainer, those corresponding to the actual stellar mass. This conclu- 
sion is consistent with our previous discussion and was obtained in 
an independent way. 

The best fit also indicates that V 393 Sco is found after a burst 
of mass transfer, the gainer having received 4 M© in a rapid burst 
lasting 400.000 years. The donor is an inflated (Rc = 9.6 Rq) and 
evolved 2 Mq star with its core consisting of 93% of helium. Ac- 
cording to the best model the system has now an age of 7.00 x 10^ 
yr and M = 9.47 x 10"^ M© yr~^ (Fig. 5). 

We notice that the mass transfer rate derived from model fit- 
ting is much smaller than those derived in Section 4.3. This could 
indicate that: (i) M values are still not well reproduced by the mod- 
els or (ii) basic assumptions of Section 4.1 are invalid, like accre- 
tion powered disc luminosity and the hypothesis of a steady state 
accretion disc. The fact that we obtain different values of M us- 
ing the accretion-powered hypothesis and that we do not observe 
orbital period changes argue for the second alternative. If viscous 
dissipation in an accretion disc is not the source of disc luminosity 
then other mechanisms should be invoked, like shocks produced in 
impact regions or photoionization of circumstellar gas and recom- 
bination of high energy photons produced in hot regions. If the disc 
is considered an extension of the gainer photosphere, it should fit its 
temperature at the inner region but to show cooler temperatures at 
the outer regions, as indicated by Eq. 2. In this case we expect the 
spectral energy distribution be characterized by thermal radiation 
from an optically thick disc in thermal contact with the hotter star. 
In this view, it is the gainer that maintains the hot disc, not viscous 
dissipation. 

The low M model imposes certain problems for the accretion 
powered disc, but explains the absence of measured orbital period 
changes. On the other hand, it is not impossible to assume that low 
M values maintain the long cycle. If we assume AI = 9 x W^'^ 
Mq yr^^ then 6.2 x lO^'* Mq are ejected in steady state in a long 
cycle. The low M solution implies that V 393 Sco is in a relatively 
long-lasting evolutionary stage compared with the previous rapid 
mass transfer rate stage. This implies that a relatively large num- 
ber of systems should be found showing the DPV phenomenon, as 
actually observed. 

If we require the best fit model to match the high M value 
then we should look at the end of the rapid mass transfer burst 
(younger system with age 6.8 x lO'^ years) or at the beginning 
of the following mass transfer event (older system with age 7.55 x 
10*^ years. Fig. 5). In both cases we can get M values of the order of 
10"'' M0 yr"^. However, these models have larger values than 
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Figure 3. Evolutionary tracks for single stars with solar metallicity (Claret 
2004) and the position for the stellar components of V 393 Sco. Derived 
masses and log g values are given between parenthesis. The evolutionary 
tracks are labeled with initial masses. 




1,5 I I 

4,8 47 4,6 4,5 4,4 4,3 4,2 4,1 40 3,9 3,8 3,7 

log L„(K) 

Figure 4. Evolutionary tracks for the binary star model from Rensbergen 
et al. (2008) that best fit the data. Donor (right track) and gainer (left track) 
evolutionary paths are shown, along with the observations for V 393 Scorpii 
(with error bars and connecting line). The best fit is reached at the time 
con'esponding to the model indicated by large crosses, that is characterized 
in Table 3. The mismatch for the primaiy is discussed in the text. Stellar 
sizes are proportional to the circle diameters. 

the preferred model by a factor of 10. In both cases, our conclusion 
that a massive structure surrounds the gainer remains unaltered. 

4.5 Double Periodic Variables in the Milky Way 

We performed a search for Galactic DPVs in the ASAS catalog 
using the software described in Section 3. Light curves were in- 
spected for the presence of two periodicities. We found the 1 1 sys- 
tems shown in Table 4 (see also Mennickent et al. 2009 and Michal- 
ska et al. 2010 for previous reports on this investigation). Two of 
our detected DPVs were previously cataloged as Algol type vari- 
ables with additional long-term variability, viz. AU Mon and V393 
Sco. The lack of dedicated stellar variability surveys in our Galaxy 
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Figure 5. M2 (upper curve) and M2 for the best evolutionary model. The 
vertical dashed line indicates the position for the best model, is shown 
in the inset graph. 
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Figure 6. The Double Periodic Variables listed in Table 4 (solid dots) and 
the hot emission-line binaries with cyclic long-term brightness changes 
listed in Table 5 (open squares). The upper and lower lines represent the 
best fit for 30 SMC DPVs (Mennickent et al. 2003) and 125 LMC DPVs 
(Poleski et al. 2010). The position of V 393 Scorpii is indicated by dotted 
lines. 



probably explains why the long photometric variability of many 
of these DPVs remained undetected for a long time. In Table 5 
we also show those Galactic hot emission-line binaries with cyclic 
long-term brightness changes reported by Desmet et al. (2010, their 
Table 10). 

The hallmark of the DPV phenomenon, i.e. the correlation 
between the orbital and the long cycle is shown in Fig. 6. V 393 
Sco fits very well the general tendency. However, two of the 
objects Usted by Desmet et al. (2010), namely /3 Lyr and RXCas 
clearly deviate from it. Excluding these two objects, the best linear 
fit for the remaining 13 systems is: 

Plong = 32.7(9)P„ (12) 

This relation is comparable to previous relationships found for 



LMC and SMC DPVs, with respective coefficients 33.1 and 32.4 
(Desmet et al. 201 1 and Mennickent et al. 2005). 

According to this relationship, the predicted long period for 
/? Lyr is about 423 days, 1.50 times longer than the reported 282 
days, and for RXCas is 1057 days, i.e. 2.05 times longer than re- 
ported. The commensurability of predicted and observed periods is 
notable. In order to investigate if the expected 423 d periodicity is 
present in the /? Lyr photometry, we removed the orbital variabil- 
ity from the V-band light curve kindly provided by Dr. Harmanec 
and searched for additional frequencies. The residual periodogram 
shows the already known periods of 282 and 356 d but only a very 
small peak in 426 days. Residuals folded with this period show a 
very noisy diagram. 

It is possible that a different evolutionary stage explains the 
position of the outliers /3 Lyr and RX Cas in Fig. 6. These objects 
could still be inside the burst of mass transfer illustrated in Fig. 5, 
while DPVs should already have passed this stage. In favor of this 
interpretation is the fact that orbital period changes have been de- 
tected in j3 Lyr and RX Cas and interpreted as due to very large 
mass transfer rates (~ 10~^ M© yr~^). However, orbital period 
changes have not been reported in DPVs. 



5 CONCLUSIONS 

In this paper we have modeled the orbital light curve of the 
intermediate-mass interacting binary V 393 Sco to obtain stellar 
and system parameters. We also disentangled the long-term light 
curve at optical and infrared photometric bands. We have found 
insights on the system evolutionary stage and long cycle nature. 
The main results of our research are: 

• The orbital period change, if present, is shorter than 0.5 sec 
per year. 

• The long-term light curves are characterized by a smooth os- 
cillation in a time scale of 253 days with larger amplitudes in redder 
bandpasses. 

• The best fit to the orbital light curve requires a non-stellar 
component that was modeled with an optically thick disc model. 
The disc radius is about half of the Roche lobe radius of the gainer 
and two bright spots are required to fit the observations. 

• We found the stellar and system parameters that best match 
the observations, which are given in Table 2 along with parameters 
for the disc and the bright spots. 

• The stability of the orbital light curve suggests that the stellar 
+ disc configuration remains stable during the long cycle. Variabil- 
ity of the optically thick disc is not the main source for long cycle. 

• Therefore and in order to fit the redder color at long maximum, 
we argue that the long cycle is produced by free-free emission in 
a variable structure, probably visible perpendicular to the orbital 
plane, something reminiscent of the jets found in /3 Lyr (Harmanec 
et al. 1996). The suggestion of equatorial mass loss as the cause of 
the long cycle by M 10 is probably biased by detection of equatorial 
outflows not necessarily related to the long cycle. 

• A comparison with published evolutionary tracks provides an 
estimate for the age of the system, namely log t = 7 x 10*^ yr. 
We find the system after a mass exchange episode, where 4 M© 
were transferred from the donor to the gainer in a period of 400.000 
years. 

• The evolutionary model with initial stellar masses of 6 and 
3.6 M© reproduces relatively well the present donor parameters 
and orbital period, but overestimates the gainer temperature and 
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Table 4. Double-periodic vaiiables in the Milky Way detected from the analysis of the ASAS database. Remarks are from SIMBAD jhttp://simbad.u-| 
^ strasbg.fr/simbad/ 1 except for V 393 Sco that are from this paper. 



Names RA(2000) DEC(2000) P^rb Plong TYCHO-2 Remarks 

HD TYCHO GCVS (d) (d) Bt Vt 



50526 


161-1014-1 




06:54:02.04 


-1-06:48:48.5 


6.7015 


191.7 


8.342 


8.244 


B9 


50846 


54801-1012-1 


AUMon 


06:54:54.71 


-01:22:32.9 


11.1132 


419.0 


8.453 


8.432 


Eel B3V-HF811Ie 




5978-472-1 




07:26:41.41 


-22:08:53.7 


8.2962 


311.0 


10.641 


10.585 






5985-958-1 




07:44:15.30 


-17:58:45.6 


7.4062 


229.5 


10.713 


10.641 






8175-333-1 


DQVel 


09:30:34.22 


-50:11:54.0 


6.0833 


188.0 


11.473 


10.977 


Eel 


90834 






10:27:41.61 


-59:17:04.9 


6.8148 


230.7 


9.591 


9.500 


B5Ill/lVe 


135938 


8695-2281-1 




15:20:08.44 


-53:45:46.5 


6.6477 


231.2 


9.470 


9.267 


B5/B61Vp 




8708-412-1 


GKNor 


15:34:50.92 


-58:23:59.1 


6.5397 


221.0 


11.646 


11.275 


Eel 


328568 


8325-3366-1 


LPAra 


16:40:01.78 


-46:39:34.9 


8.5331 


274.0 


10.541 


10.243 


EclB8 


161741 


7385-1101-1 


V393 Sco 


17:48:47.60 


-35:03:25.6 


7.7126 


251.7 


7.748 


7.609 


Eel B4V/A7III 


170582 


5703-2382-1 




18:30:47.53 


-14:47:27.8 


16.872 


536.0 


10.150 


9.711 


A3-A9 



Table 5. Galactic hot emission-line binaries with cyclic long-term brightness changes (Desmet et al. 2010). P^^ is the predicted long period according to 
Eq. 12. References for the spectral types are given. 





Names 




RA(2000) 


DEC(2000) 


-Porb 




Pth 


Remarks 


HD 


TYCHO 


GCVS 






(d) 


(d) 


(d) 






4313-258-1 


RXCas 


03:07:45.75 


-1-67:34:38.6 


32.312 


516 


1057 


A5III-I-G3111: Strupat 1987 


174237 


3918-1829-1 


CXDra 


18:46:43.09 


-1-52:59:16.7 


6.696 


130-180 


219 


B2.5V+F5m; Simon 1996 


174638 


2642-2929-1 


/3Lyr 


18:50:04.80 


-1-33:21:45.6 


12.94 


282.4 


423 


B8eplI-HB6.5; Budding et al. 2004 


216200 


3223-3619-1 


V360 Lac 


22:50:21.77 


-1-41:57:12.2 


10.085 


322.2 


330 


Be-l-F; Linnell et al. 2006 



luminosity, a fact that could be ascribed to the optically thick disc 
not considered in the evolutionary tracks. In order to explain these 
features, we argue that part of the mass, maybe up to 2 M© has not 
been accreted by the gainer, but remains in the massive optically- 
thick disc. 

• We find a discrepancy between values of Al derived from 
theoretical model fitting of observationally-derived parameters and 
those derived from accretion-theory analysis of observationally- 
derived parameters, the former being smaller by 2 orders of mag- 
nitude. If the optically thick disc is a representation of a mas- 
sive disc-like pseudo-photosphere, then its luminosity could not be 
accretion-driven. The constancy of the orbital period supports the 
view that the disk luminosity is not driven by viscosity, but proba- 
bly by reprocessed stellar radiation. 

• We present the results of our search for Galactic DPVs and 
make a comparison with hot emission-line binaries with cyclic 
long-term brightness changes. 13 Galactic DPVs show a similar 
correlation between Po and Piong to that observed in LMC and 
SMC DPVs. The systems /3 Lyr and RX Gas deviate from this ten- 
dency. These systems could be in an earlier evolutionary stage com- 
pared with DPVs. 
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